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ABSTRACT 

We follow the structural evolution of star-forming galaxies (SFGs) like the Milky Way by selecting progenitors to 
z ~ 1.3 based on the stellar mass growth inferred from the evolution of the star-forming sequence. We select our 
sample from the 3D -HST survey, which utilizes spectroscopy from the HST/ WFC3 G141 near-IR grism and enables 
precise redshift measurements for our sample of SFGs. Structural properties are obtained from Sersic profile fits to 
CANDELS WFC3 imaging. The progenitors of z = 0 SFGs with stellar mass M = 10 10,5 M© are typically half 
as massive at z ~ 1. This late-time stellar mass growth is consistent with recent studies that employ abundance 
matching techniques. The descendant SFGs at z ~ 0 have grown in half-light radius by a factor of ~1 .4 since z ~ 1 . 

The half-light radius grows with stellar mass as r e a M° 29 . While most of the stellar mass is clearly assembling at 
large radii, the mass surface density profiles reveal ongoing mass growth also in the central regions where bulges 
and pseudobulges are common features in present day late-type galaxies. Some portion of this growth in the central 
regions is due to star formation as recent observations of Ha maps for SFGs at z ~ 1 are found to be extended but 
centrally peaked. Connecting our lookback study with galactic archeology, we find the stellar mass surface density 
at R = 8 kpc to have increased by a factor of ~2 since z ~ 1, in good agreement with measurements derived for 
the solar neighborhood of the Milky Way. 
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1. INTRODUCTION 

The structural formation of late-type star-forming galaxies 
(SFGs) like the Milky Way has been studied through several 
complementary approaches. Detailed studies of the ages, metal- 
licities, and kinematics of stellar populations within the Milky 
Way itself afford a unique vantage point for viewing the assem- 
bly of such late-type galaxies — although such an approach is 
limited in sample size and completeness (see, e.g., Rix & Bovy 
2013). On the theoretical front, the formation of Milky -Way-like 
SFGs has posed a challenge for simulations that aim to produce 
a structurally realistic analog at z ~ 0 with a disk and an embed- 
ded bulge or pseudobulge (Scannapieco et al. 2012). This largely 
reflects the difficulty of modeling the “sub-grid” physical pro- 
cesses that impact the baryons (although see, e.g., Guedes et al. 
2011; Stinson et al. 2013). Lookback studies provide another 
window into the formation of these galaxies. Though several 
works have studied the structural properties of late-type galaxies 
in the distant universe, including the size-mass relation to z ~ 1 
and beyond (e.g., Barden et al. 2005; Trujillo et al. 2006), the 
connection to the structural evolution of a typical disk galaxy, as 
it grows in mass, remains largely unexplored (for comparisons 
to models, see Somerville et al. 2008; Dutton et al. 2011). 

In this paper, we systematically select progenitors of Milky- 
Way-like SFGs using the stellar mass growth inferred from 


* Based on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute. STScI is operated by the 
Association of Universities for Research in Astronomy, Inc. under NASA 
contract NAS 5-26555. 


the evolution of the star-forming sequence and analyze their 
structural evolution. We use HST /WFC3 near-IR imaging that 
probes redder rest- frame wavelengths than most previous works 
and is therefore less sensitive to light from young stars. Our 
selection provides a quantitative view for stellar mass build-up 
in different radial regimes. The method for computing the mass 
growth of SFGs has been discussed in detail by Leitner (2012) 
and complements other studies that connect progenitors and 
descendants using number densities (van Dokkum et al. 2010; 
Patel et al. 2013). Given the distinct formation history of SFGs 
from quiescent galaxies (QGs), we use the former method here 
as it directly traces the star-forming progenitors of galaxies like 
the Milky Way. In van Dokkum et al. (2013), we use the number 
density approach to study progenitors of galaxies of all types 
with the mass of the Milky Way. 

We assume a cosmology with Hq = 70 km s _1 Mpc -1 , 
Q M = 0.30, and Qa = 0.70. Stellar masses are based on a 
Chabrier initial mass function (Chabrier 2003). Magnitudes are 
in the AB system. 

2. DATA AND ANALYSIS 

We employ data from the 3D-HST survey (v2.\) to carry 
out our analysis. The observations and data reduction proce- 
dures are explained in detail in Brammer et al. (2012) and R. E. 
Skelton et al. (in preparation). The HST/ WFC3 G141 near-IR 
grism observations are the centerpiece of 3D -HST and cover the 
CANDELS fields. In this work, we use the three fields with cur- 
rently available WFC3 based structural parameters from the 


1 


The Astrophysical Journal, 778:115 (7pp), 2013 December 1 


Patel et al. 



Figure 1 . (a) Evolution of SFGs in the SSFR-mass plane with final mass at z ~ 0 of M = 10 105 Mq as computed by Leitner (2012) using the SSFR-mass relations 
at different redshifts (dashed lines) from Karim et al. (2011). At a given redshift, new stellar mass is added according to the position of the SFG on the star-forming 
sequence and mass loss from stellar evolution is also accounted for. The slight offset between the extrapolated z = 0 relation from that of Salim et al. (2007) does not 
impact our results (see text). The Milky Way (gray circle) lies within the observational scatter of the z = 0 relation, (b) SFH of the SFGs tracked in panel (a). The 
data points with error bars represent IR+UV SFRs for a subset of our sample with deep MIPS: there is good general agreement with the radio based SFRs from Karim 
et al. (2011). (c) Mass growth history of the SFGs tracked in panel (a). We select galaxies at the indicated masses and redshifts in this work (i.e., colored points on 
the black curve). Since z = 1, the SFGs grow in stellar mass by a factor of ~2. Mass growth histories derived from alternative SFR measurements give very similar 
results (e.g., green curve, which uses far-infrared SFRs from Oliver et al. 2010). For comparison, the SFH and mass growth for galaxies with the same final mass from 
the abundance matching work of Behroozi et al. (2013) is shown by the dashed red line. There is good agreement between the two different methods. 

(A color version of this figure is available in the online journal.) 


literature: COSMOS, GOODS-S, and UDS. Grism redshifts 
were measured using a modified version of EAZY (Brammer 
et al. 2008). The procedure fits template spectral energy distribu- 
tions to both the photometry and grism spectroscopy, enabling 
precise redshift measurements (Brammer et al. 2012) for our 
sample of SFGs at high redshift. Objects selected in this work 
had full spectral coverage and less than <50% integrated spec- 
tral contamination. 

Stellar masses were measured with FAST (Kriek et al. 2009). 
We determine the stellar mass limit at our highest redshifts 
(z = 1.4) to be M ~ 10 9 - 8 Mq (see Figure 1(c)) using a similar 
technique to that of Marchesini et al. (2009). This conservative 
limit accounts for galaxies with high M/L such as reddened 
SFGs. Rest-frame U—V and V — J colors were measured with 
EAZY in order to distinguish SFGs from QGs (Section 3.2). 

Structural properties were obtained from van der Wei et al. 
(2012) who used GAFFIT (Peng et al. 2010) to fit Sersic profiles 
to the CANDEFS HST/ WFC3 J 125 and H m imaging (Grogin 
et al. 2011; Koekemoer et al. 2011). Only objects with no flags 
were selected. The effective radii reported here are circularized. 
The J \ 25 imaging is used at 0.25 < z < 1 while the imaging 

is used at 1 < z < 1.4, with our median galaxy sampling rest- 
frame ~7100 A. The point-spread function (PSF) FWHM/2 of 
the //160 imaging is MX'09, corresponding to a physical radius 
at our highest redshifts of ~0.8 kpc. 

We employ a Sloan Digital Sky Survey (SDSS) reference 
sample at z = 0.05 using stellar masses and star formation rates 
(SFRs) from the MPA/JHU catalogs (Brinchmann et al. 2004) 
and i- band Sersic fits from the NYU-VAGC (Blanton et al. 2003, 
2005). 

3. SEFECTION 

3.1. Stellar Mass Growth History from 
an Evolving Star-forming Sequence 

In this work, we trace the formation of SFGs with a final mass 
of M — 10 10,5 Mq, which is slightly below (~0.2 dex) that of 
the Milky Way, so that we can compare various properties to 
galaxies of a similar final mass from other works that employ 


different progenitor-descendant linking methods (e.g., Behroozi 
et al. 2013). We note however that we arrive at similar qualitative 
conclusions when tracing progenitors of Milky-Way-mass SFGs 
(i.e., M z= 0 ~ 10 10 - 7 M q ). 

Selecting progenitors of z ~ 0 SFGs requires knowledge of 
their mass growth history so that one can select galaxies of the 
proper progenitor mass at a given redshift. For galaxies that 
assemble most of their stars from in situ star formation, one 
can infer the mass growth from the evolution in the SFR-mass 
relation (or SSFR-mass). The Milky Way, with a stellar mass 
of M - 5 x 10 10 Mq (Hammer et al. 2007) and a SFR of 
1.9 ± 0.4 Mq yr -1 (Chomiuk & Povich 2011), falls within 
the observational scatter of the SSFR-mass relation of Salim 
et al. (2007, their Equation (11), corrected for evolution to 
z = 0 assuming Equation (1) here; see Figure 1(a)). This star- 
forming sequence has been extensively studied to high redshifts 
(e.g., Noeske et al. 2007; Karim et al. 2011; Whitaker et al. 
2012; Fumagalli et al. 2012). The method for determining the 
mass growth from this relation is simple: for a given redshift 
interval, new stellar mass is added to the existing mass based 
on the location of a galaxy in the SSFR-mass plane and mass 
loss from stellar evolution is accounted for from simple stellar 
population models. The assumption that most nearby SFGs were 
star forming at high redshift is supported by the small scatter 
about the SSFR-mass relation as well as by the vastly different 
structural properties between SFGs and QGs (e.g., Franx et al. 
2008). 

In this work, we use the mass growth computed by Feitner 
(2012), who derived it from the SSFR-mass relations of Karim 
et al. (2011). Figure 1(a) shows SSFR-mass relations at different 
redshifts (dashed lines) measured by Karim et al. (2011) from 
1.4 GHz radio stacking. They employ the parameterization 

SSFR oc M^(l + z) n (1) 

where f> = —0.35 and n — 3.45. Though the z = 0 
relation based on Equation (1) is an extrapolation of the Karim 
et al. (2011) data, its slight difference from lower redshift 
observations (e.g., Salim et al. 2007) does not significantly 
impact the mass growth at these late times. The track shown 
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Figure 2. ((a)-(e)) Rest-frame U—V vs. V—J at different redshifts, with galaxies selected to have the progenitor mass (±0.1 dex) from Figure 1(c). The Williams 
et al. (2009) boundary separating QGs from SFGs (bottom right) is used to select the latter. The gray curve indicates the evolution in UVJ colors for the SFH shown in 
Figure 1(b) with the darker segment representing the colors spanning the given redshift interval. The SFGs are color-coded according to their Sersic indices. Sample 
sizes are indicated in the bottom right, (f) WFC3 G141 median stacked spectra of SFGs for redshift bins with observable Ha. 

(A color version of this figure is available in the online journal.) 


in Figure 1(a) (black curve) indicates the trajectory in the 
SSFR-mass plane for the SFGs selected for study in this 
work with final mass at z = 0 of M = 10 10,5 M o . The 
corresponding star formation history (SFH) is shown in panel 
(b) and the mass growth with redshift in panel (c). The points 
with error bars (including a 30% systematic uncertainty) in 
panel (b) indicate the median SFR (IR±UV, computed similarly 
to Franx et al. 2008) of our sample in the GOODS fields (i.e., 
deep MIPS) for our LV/-selected (see below) SFGs at the 
corresponding redshifts and masses from panel (c). While SFR 
estimates can vary due to systematics between different tracers 
(e.g., Wuyts et al. 2011; Leitner 2012) there is good general 
agreement between our SFR measurements and those of Karim 
et al. (2011). The green curve in panel (c) indicates the mass 
growth computed by Leitner (2012) when employing the far- 
infrared derived SFR data of Oliver et al. (2010; f> = —0.15, 
n = 3.36). The small difference (<0.05 dex) compared to the 
adopted mass growth from the Karim et al. (2011) SFR data 
(black curve) indicates that variations in the sample selection 
due to uncertainties in the mass growth from different SFR 
data are small and therefore does not impact our qualitative 
conclusions. The dashed red curves in panels (b) and (c) show 
the corresponding SFH and mass growth for the same final 
mass galaxy from the abundance matching work of Behroozi 


et al. (2013). There is good agreement between the two different 
methods. 

The results above imply significant stellar mass growth in our 
SFGs below z < 1 (factor of ~2.2). Late mass growth at z < 1 
is also found by other works (e.g., Yang et al. 2012; Moster 
et al. 2013). SFHs of the Milky Way disk derived from stellar 
properties also suggest significant late-time assembly (e.g., 
Rocha-Pinto et al. 2000; Aumer & Binney 2009). Finally, disks 
are more commonly formed in hydrodynamical simulations that 
favor late-time assembly (Scannapieco et al. 2012). 

3.2. Selection of Star-forming Progenitors to z ~ 1.3 

Given the observational scatter about the SSFR-mass rela- 
tion, in practice we select all SFGs in a given redshift and mass 
bin. Figure 2 shows rest- frame U—V versus V — J diagrams for 
galaxies in narrow mass bins (±0.1 dex) centered on the pro- 
genitor mass at a given redshift. We use this UVJ selection to 
identify SFGs, which occupy the bottom right portion of these 
diagrams (Williams et al. 2009; Patel et al. 2012). The gray 
curve in the figure indicates the evolution in UVJ colors for the 
dust-free SFH shown in Figure 1 (b) with the darker segment rep- 
resenting the colors spanning the given redshift interval. These 
UVJ colors are extended redward depending on the amount 
of reddening from dust (see, e.g., the reddening vector in the 
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Figure 3. Random SDSS i - band and HST/ WFC3 postage stamps for progenitors of SFGs with a final mass at z = 0 of M = 10 105 Af 0 . Each stamp is ~30 kpc on a 
side. The median stellar mass decreases to high redshift according to Figure 1(c) where the sizes appear smaller. 

(A color version of this figure is available in the online journal.) 


figure). The descendants, which are more massive, become red- 
der toward low redshift likely due to their aging stellar popu- 
lations. We note that the Williams et al. (2009) boundary that 
is used to distinguish SFGs from QGs shifts to redder U—V 
colors toward low redshift and therefore accounts for the gen- 
eral decline in galaxy SFHs. In the SDSS sample, galaxies with 
SSFR >10 11 yr 1 were selected as SFGs. Figure 3 shows ran- 
dom SDSS i- band and //5T/WFC3 postage stamps for progen- 
itor SFGs at different redshifts. 

4. RESULTS 
4.1. Structural Evolution 

We examine the evolution in various structural parameters for 
our sample of progenitor SFGs in Figure 4. Panel (a) shows that 
most of the progenitors at high redshift have close to exponential 
profiles while those at lower redshift have slightly higher Sersic 
indices. The median axis ratio in panel (b) remains relatively 
constant at b/a ~ 0.60, a low value which for a population 
of randomly inclined disks modeled as oblate spheroids would 
imply an intrinsic axis ratio of ~0.3. The low Sersic indices and 
low axis ratios for the progenitor SFGs at higher redshifts is 
suggestive of disks. 


Half-light radii provide a first order view into the distribution 
of stellar mass for galaxies in our sample. Figure 4(c) shows the 
evolution of the median half-light radius with redshift. Though 
we follow the median of a given property, it is important to note 
that at a given redshift and stellar mass, SFGs display a diversity 
of property values. The constant scatter with redshift about 
the size-mass relation (A. van der Wei et al., in preparation), 
however, suggests that the increasing intercept of this relation 
toward low redshift is driven by evolution in the SFG population 
as a whole. The black dashed line represents a fit to the data of 
the form 

r e = P(l+z) a (2) 

where /3 = 3.5 zb 0.3 kpc and a = —0.48 ± 0.15. We note 
that all of the SDSS data points in Figure 4 include a 10% 
error added in quadrature to account for systematics between 
different measurement methods (e.g., Figure Al in Guo et al. 
2009). Since z = 1, the median half-light radius for these SFGs 
has grown by a factor of ~1.4, indicative of some amount of 
inside-out growth. 

Figure 4(d) shows the evolution of the progenitors in the 
size-mass plane. For reference, the SDSS z ~ 0 size-mass 
relation for late-type galaxies from Shen et al. (2003) is indicated 
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Figure 4. (a) Median Sersic index vs. redshift for progenitor SFGs. (b) Median axis ratio vs. redshift. (c) Median half-light radius vs. redshift. The error bars represent 
the bootstrapped uncertainty on the median. Fits of the form (1 + z)“ are indicated for the Sersic indices and half-light radii, both of which increase toward low redshift. 
(d) Evolution of r e in the size-mass plane. The dotted blue line is the z ~ 0 SDSS size-mass relation (z-band) for late-type galaxies from Shen et al. (2003). The 
dashed line is a fit to the data of the form r e oc M a , where a = 0.29 ± 0.08. Since z = 1, the half-light radii of the SFGs have grown by a factor of ~1.4. 

(A color version of this figure is available in the online journal.) 


by the dotted blue line. The dashed line represents a fit to our 
sample of the form: 

r e oc M a (3) 

where a = 0.29 =b 0.08. Interestingly, the value of a measured 
here for SFGs is much smaller than that for QGs (o:qg ~ 2; van 
Dokkum et al. 2010; Patel et al. 2013). SFGs generally appear to 
evolve near the local scaling relation below z < 1.3, consistent 
with slow evolution in the intercept of the size-mass relation for 
SFGs (A. van der Wei et al., in preparation). Finally, we note that 
none of the results presented in Figure 4 change significantly 
when using the mass growth derived from the Oliver et al. (2010) 
SFRs (green curve in Figure 1(c)). 

4.2. Stellar Mass Growth in the Central 
and Outer Regions at z < 1-3 

Stellar mass surface density profiles provide a more detailed 
look at the distribution of mass within galaxies. Figure 5(a) 
shows the median combined mass surface density profiles for 
our sample of SFG progenitors. These profiles were constructed 


in a similar manner to that of Patel et al. (2013), where the 
best-fit single component Sersic profiles were stacked to create 
the median light profile. A single component Sersic profile is 
generally found to be a good representation of the light profile 
for individual galaxies at high redshifts (Szomoru et al. 2012; 
Mosleh et al. 2013). Prior to stacking, each light profile was 
converted into a mass profile by normalizing the light within 
R = 20 kpc to the total stellar mass of each galaxy. While a 
more thorough analysis would be required to properly account 
for M/L gradients (see, e.g., Szomoru et al. 2013), the present 
analysis provides a first order view into the evolution in the radial 
distribution of stellar mass. Below, however, we also examine 
how different M/L values in the central and outer radial regimes 
impact the measured stellar mass growth. In addition, we will 
explore the impact of M/L gradients for a broader sample in 
followup work. Figure 5(b) shows the ratio of the mass profiles 
at higher redshifts to the SDSS mass profile. Below z < 1.3, 
stellar mass is continually built up at all radii. 

Figure 5(c) shows the growth in the projected mass inside 
(red) and outside (blue) of R = 2 kpc (i.e., ~ r e at z ~ 
1.3). Clearly more mass has assembled at larger radii since 
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Figure 5. (a) Median combined mass surface density profiles. The dotted portion of each curve indicates where the uncertainty becomes larger than >20%. The shaded 
gray region indicates the maximum WFC3 PSF FWHM/2. The solar symbol indicates the mass surface density in the solar neighborhood as computed by Bovy et al. 
(2012) but scaled down by 0.2 dex to account for the difference in mass between the Milky Way and our z ~ 0 descendants, (b) The ratio of the mass profiles at higher 
redshifts to the SDSS mass profile. Shaded region indicates the uncertainty. The factor of 2.4 ± 0.3 mass growth from z ~ 0.9 to z ~ 0 at Ro = 8 kpc agrees well with 
measurements in the solar neighborhood from Aumer & Binney (2009). (c) Median stellar mass growth for the central (R < 2 kpc; red) and outer regions (R > 2 kpc; 
blue). The total stellar mass growth is shown in black. While most of the new mass growth since z ~ 1 has taken place in the outer regions, mass has also been added 
in the central regions of these SFGs where bulges and pseudobulges are common features. Some portion of the growth in the central regions is due to star formation as 
Ha maps for SFGs at z ~ 1 are centrally peaked (Nelson et al. 2013). Dashed lines with open circles indicate the mass growth when accounting for M/L differences 
in the two radial regimes (see text). 

(A color version of this figure is available in the online journal.) 


z ~ 1.3, leading to larger r e toward low redshift. However, the 
amount of stellar mass in the central regions, where bulges and 
pseudobulges are common structural features in nearby late-type 
galaxies (Weinzirl et al. 2009), has also increased. Some portion 
of the growth in the central regions is due to star formation as Ha 
maps for SFGs at z ~ 1 are centrally peaked (Nelson et al. 2013). 
This is qualitatively consistent with recent observations of stars 
in the Milky Way bulge that display a wide range of metallicities 
and ages, implying an extended formation history (Bensby et al. 
2013). The ongoing mass assembly in the central regions may 
point to secular processes (e.g., Kormendy & Kennicutt 2004) 
as an important channel for bulge growth for SFGs in the stellar 
mass regime studied here. 

We estimate the impact on the stellar mass growth in the two 
different radial regimes due to variation in M/L ratios. These 
M/L ratios can be derived from the correlation with galaxy 
colors (e.g., Bell & de Jong 2001). For each galaxy in our 
sample, we use the PSF deconvolved best-fit GALFIT models 
in the Jus and bands to estimate the observed Jus — #i6o 
colors inside and outside of R = 2 kpc (R = 10 kpc was used for 
the upper bounds on the outer radial regime). The median color 
difference between the different radial regimes (outer minus 
central) is, from low (0.25 < z < 0.5) to high redshift (1.2 < 
z < 1.4), A (Jus ~ H m ) = -0.09,-0.11,-0.13,-0.10, 
and —0.09 mag. The outer regions are therefore bluer at all 
redshifts. The A(Jus — #i6o) values were then converted into 
A(log M/ L x ) values using the slope of the correlation between 
log M/L x and Jus — Hm, where X represents the observed 
J us (for the sample at z < 1) and Hm (z > 1) bands. In 
deriving these M/L - color correlations at different redshifts, we 
employed a broad sample of galaxies, including both QGs and 
SFGs above a mass of M > 10 9 5 M© so that a wide array of 
stellar populations were sampled. From low to high redshift, 


we find A(log M / L x )/A( Jus - H m ) = 0.44, 1.33, 1.74, 1.33, 
and 1.45 dex mag -1 . The scatter in lo gM/L x about these 
correlations is 0.15-0.19 dex. While some of this scatter is 
caused by working in the observed frame with a sample spanning 
a range of redshifts, it is similar to what is found in other recent 
works (e.g., 0. 13-0.28 dex in Szomoru et al. 2013). We multiply 
our median light profile at a given redshift by the step function 
implied by A(lo gM/L x ) (re-normalizing to the median mass 
of the sample) and integrate as before to determine the mass in 
the central and outer regions (open circles with dashed lines in 
Figure 5(c)). For the SDSS light profile, we apply the i- band 
M/L gradient computed by Tortora et al. (2011, their Figure 3). 
We adopt the average value between their early-type galaxy and 
late-type galaxy samples (A(log M/L *)/A(log R/R e ) ~ —0.1) 
since these two classes were divided by, among other properties, 
their Sersic indices with a cut at n = 2.5, which is close to the 
median value for our SDSS sample. The overall conclusion, 
compared to the case where M/L is assumed to be constant 
with radius, is unchanged, as both radial regimes undergo mass 
growth since z ~ 1 but with the outer regions building up 
relatively more mass. 

5. DISCUSSION 

5.1. A Comparison with the Stellar Mass Growth 
in the Solar Neighborhood 

Given the significant mass assembly found at large radii, 
we compare our results in such a region of our galaxy that 
has been well documented, the solar neighborhood. Though we 
caution that the Milky Way is just one such SFG, and may not 
be one that is an archetypal late-type (Hammer et al. 2007), 
this analysis nevertheless provides an intriguing comparison 
between our lookback study and galactic archeology. The solar 
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symbol in Figure 5(a) indicates the stellar mass surface density 
at the solar radius ( Ro = 8.0 kpc; Vanhollebeke et al. 2009) 
for the Milky Way from Bovy et al. (2012) but scaled down by 
0.2 dex to account for the difference in mass with our descendant 
SFGs at z ~ 0. The systematic uncertainty in this correction 
may account for the slight offset from the SDSS mass profile, 
though the rough agreement is still remarkable. Aumer & Binney 
(2009) estimate a mean formation time for stars in the solar 
neighborhood that corresponds to z ~ 0.9. Assuming this mean 
is close to the median formation redshift, therefore implying a 
factor of ~2 growth below z < 0.9, this estimate is close to the 
factor of ~2.4 zb 0.3 growth in our mass surface density from 
z ~ 0.9 (green line) to the SDSS mass profile at R = 8 kpc. 

5.2. Caveats 

Two caveats to the analysis presented here warrant some 
consideration. First, the contribution of stars formed ex situ 
to the stellar mass growth is an uncertain quantity, though there 
is evidence that it is minimal. Both Behroozi et al. (2013) and 
Moster et al. (2013) find that for local halos of mass M h ~ 
10 12 Mq (hosting centrals of stellar mass M = 10 10,5 AT 0 ), little 
stellar mass was accreted. At the lowest redshifts, Behroozi et al. 
(2013) find a growing contribution, but this is likely driven by the 
dominance of QGs in their halos at low redshift, which primarily 
grow from mergers. Major mergers are likely rare in our sample 
at z < 1 as such events tend to destroy disks, contrary to what 
is observed at low redshift (Figure 3). 

Second, the QG fraction increases toward low redshift and as 
a result not all SFGs at z ~ 1.3 will be progenitors of SFGs at 
low redshift as some will have quenched. This would impact our 
results if there is a relation between the structure of high redshift 
SFG progenitor candidates and low redshift QG descendants, 
as the former would bias our measurements for the median 
structural property. 

5.3. Conclusions 

In this paper, we have used the evolution in the SSFR-mass 
relation to determine the expected stellar mass growth (e.g., 
Leitner 2012) for progenitors of SFGs like the Milky Way. 
We used 3D -HST redshifts and photometry to select SFGs of 
the appropriate progenitor mass at different redshifts back to 
z ~ 1 and HST CANDELS imaging to follow their structural 
evolution. As these SFGs grew in stellar mass by a factor of 
~2 since z ~ 1 , most of the new stellar mass assembled in the 
outer regions, leading to a factor of ~ 1.4 increase in the half- 
light radius. The mass surface density profiles indicate ongoing 
stellar mass assembly also in the central regions where bulges 
and pseudobulges are common features in spirals. We also found 
good agreement between the mass growth at R = 8 kpc in 
our lookback study and that for the solar neighborhood of the 
Milky Way. 

We thank Daniel Szomoru and Sean McGee for helpful 
discussions. 
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